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COMPUTATION OF TIDES AND TIDAL CURRENTS— 
UNITED STATES PRACTICE 


H.A. Einstein,* M. ASCE, and R.A. Fuchs** 


The authors of this paper have recently made a survey of past and present 
calculation methods used for the prediction of tidal stages and flows in canals 
and estuaries. The results of this study conducted under contract for the 
Committee on Tidal Hydraulics, Corps of Engineers, U.S. Army, are in part 
offered in this paper. 

As the title says, this report is limited to the U.S. Practice, as developed 
and applied before the second World War. Very little work was done in this 
field during and immediately after the war until the Committee on Tidal 
Hydraulics revived the interest, a result of this effort being the present sym- 
posium, As no new calculation methods were developed and tested which 
could be termed original during this study, the present paper will be restrict- 
ed to the description and discussion of the prominent older methods, where 
there will be some occasion to indicate where the methods may be modified 
for added accuracy or simplicity. 

Before these methods are discussed in detail, it appears to be necessary 
to introduce some of the concepts and the physical background of all such cal- 
culations. The general problem of a tidal flow with or without a constant flow 
component is mathematically an extremely complex problem and can not be 
solved rigorously in most cases. Most solutions depend on the introduction of 
various simplifying assumptions which may or may not be permissible in any 
particular application. A short resume may be given first of all the compli- 
cating aspects of the tidal flow problem in general. The calculation methods 
referred to later will be shown to describe only partially the complicated 
flows. 


The General Tidal Flow and the Tidal Basin 


Most people associate with the term “tidal flow” reversing flow caused by 
tidal action. It can easily be seen that any such flow is always affected by 
other influences, too, such as friction or the inflow of a river. It is much 
more appropriate, therefore, to include into our problem all flows which are 
to some degree affected by tidal action, which is the wave condition generated 
in large bodies of water by the variation of the attraction by the moon and the 
sun as these celestial bodies change their relative position. From this we 
may conclude that a large body of water, an ocean, must be part of the prob- 
lem since the tidal action originates in it. For the purpose of this paper, the 
flows in the large body of water in which the tidal action originates are not 
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included in the term “tidal flow.” The problem to be studied is the flow con- 
dition in smaller bodies of water which are connected with an ocean and in- 
fluenced by the water surface variations in this ocean. They will be assumed 
to be sufficiently small that any direct tidal attraction of the sun and the moon 
on its water can be neglected as insignificant. 

It will thus not be the object of this report to discuss the methods of pre- 
dicting the tidal cycles of the water surface at any point of the oceans; these 
cycles will rather be assumed to be known and will be introduced as boundary 
conditions of the flow in the estuaries, bays and canals under consideration. 
We will assume, furthermore, that the estuary under consideration is suffi- 
ciently small that its influence on the cycle of tidal water surface elevations 
in the ocean is negligible. 

In contrast to other movements of the ocean water, such as those produced 
by wind, the tides are periodic and can be predicted well in advance for a 
given location on the basis of past observations. The fluid motion is mainly 
horizontal and hence can be shown to be the same for all particles on a verti- 
cal line. The tidal elevation at a given location can be analyzed into a series 
whose periods are those of the disturbing forces and are therefore known a 
priori. The amplitudes and phases of the component tides are determined by 
harmonic analysis of actual tidal observations over sufficiently long periods. 

The larger part of the tidal variation is due to the motion of the moon be- 
cause of its proximity to the earth. The principal effects are due to varia- 
tions in the phase of the moon, distance and declination. The relative magni- 
tude of these effects varies in different parts of the world. For the precise 
prediction of the tide at any place one may need from 20 to 30 component 
tides. The five principal components are as follows: 

M, = lunar semi-diurnal with a period of 12.42 hours 

8, = solar semi-diurnal with a period of 12 hours 

N, = larger lunar elliptic semi-diurnal with a period of 12.66 hours 

K, = luni-solar diurnal with a period of 23.93 hours 

0, = lunar diurnal with a period of 25.82 hours. 

One finds that the semi-diurnal tides prevail on the shores of the Atlantic 
Ocean, the diurnal tide on the coast of China, Alaska and the Gulf of Mexico, 
and mixed tides on the Pacific Coast of North America. 

The tide at any point in the oceans may be represented in the form 


7* Not 1, COS(aot-k,) + 72 COS (2aot—ke) +.. 


where @ is the frequency of the diurnal tide = 0.702 x 10~* sec~*. On the coast 
of the Netherlands one finds, for example, that the amplitude of the M, tide is 
approximately 80 per cent of that of the total tide. Under such conditions tidal 
flow in canals can be approximated by a single harmonic component with a 
period of about 12 hours. This type of approximation used so successfully and 
extensively by the Dutch may be puuite inadequate when the tide is of the mixed 
type such as at San Francisco. (1 

While the introduction of this variable water-surface boundary condition is 
the characteristic outside influence on tidal flows, it is by no means the only 
one, In the case of canals, for instance, it is very common that such a bound- 
ary condition exists not only at one end but at both ends and that the tidal 
cycles at the two ends are different in shape, phase and amplitude. If a sys- 
tem of canals is under consideration there may be even more openings to the 
open sea, at which the tidal cycle of the ocean provides boundary conditions 
for the tidal flow in the canal system. 
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In many tidal flow problems other types of boundary conditions must be 
satisfied. It is rather common to many estuaries that a river enters the es- 
tuary at the land side. The transition from river to estuary is often so grad- 
ual that it becomes impossible to decide where one begins and the other ends. 
The safest way of handling that problem is then to assume the boundary to be 
in the river and to introduce the river flow as boundary condition. This con- 
dition represents thus a discharge, and not a stage as that at the ocean side. 
The same condition must be introduced at a dead end of the estuary—the dis- 
charge as the flow velocity at that point is zero. 

Other boundary conditions may be introduced in other cases. For instance 
a very long canal may be approximated by an infinitely long canal, which is 
mathematically much easier to solve than the canal of finite length. 

Our tidal flow problem is now defined as the flow of water in a basin which 
is small compared to the ocean, but directly connected with it, and of which 
we know some boundary conditions. This flow wants to be described by the 
water surface elevation and the flow velocity at various points of the system 
as a function of time. Such a description permits us to predict when and 
where dangerous or undesirable stages or velocities will occur in the system. 

Most canals and estuaries have a limited water depth and the flows in them 
will incur appreciable friction. In the calculation of the flow the frictional ef- 
fects can for that reason usually not be neglected. It will be seen that it is 
the turbulent friction which complicates the mathematical description of tidal 
flows more than any other single factor. 

If a river enters the estuary there will be both sweet river water and salty 
sea water in the estuary. Density currents, i.e. the relative motion in hori- 
zontal layers of water of different density may occur in this case. However, 
there is an interesting interrelation between turbulent friction and the occur- 
rence of such density currents which will permit us to separate these two ef- 
fects. Turbulent friction, such as the friction described by Manning’s equa- 
tion or similar others, is characterized by the existence of turbulence in the 
entire cross section. This turbulence has vertical components which ex- 
change fluid vertically at such a rate that the transfer of momentum between 
different layers holds the entire shear stresses in equilibrium. Such a turbu- 
lent exchange of fluid would naturally mix water layers of different density in 
a rather short time and make density currents impossible. We can conclude 
from this that density currents are probably not an important factor where 
friction is important. In those cases where the turbulent mixing is insuffi- 
cient to break up density currents, the flow patterns will be so different from 
that usually described by the common friction laws that these friction laws do 
not apply. A combination of density currents and frictional effects is thus at 
the existing stage of our knowledge impossible. 

It is a very common occurrence that sediment is moved to the tidal flow 
basin either by a river or from the open sea. Some of the most important 
engineering problems in such basins are actually caused by sediment deposi- 
tion. The question thus arises whether or not the sediment motion as such 
has an influence on the resulting tidal flows. 

For a general appraisal! of this effect it is necessary to divide the sedi- 
ment materials into a fine part, which stays in suspension until it is finally 
deposited in some dead-water area where it remains permanently, and into a 
coarser part which moves essentially along the bottom whenever the flow is 
locally sufficienily strong, and which represents a shifting deposit of loose 
material. It is clear that both types of sediment, if continuously supplied to 
the basin, will gradually fill it. This effect is sufficiently slow and gradual, 
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however, that it has no effect on the flow at a given instance, but that it will 
only change slowly the geometry of the basin. The main instantaneous influ- 
ence on the flow comes due to the fact that the bottom is covered by such sed- 
iment material and it is thus the sediment material which defines the rough- 
ness of the bottom. The fine material in the low flow areas may be observed 
to deposit as a very smooth and regular surface which results in a very low 
friction factor. Such areas will afford the flow a minimum resistance unless 
friction is caused by any foreign matter such as a protruding rock, structures 
or debris deposited at extremely large flows. The shifting coarser materials 
on the other hand are rougher in themselves. In addition they often have the 
tendency to form ripples and bars due to flow and wave action, such that they 
provide not only for a rougher, but also for a more irregular surface—both 
factors increasing the effective roughness of the bottom. A careful choice 
will thus be necessary of friction factors taking into account the distribution 
of various sediment deposits. In most basins the effective value of such fric- 
tion factors is actually possible only by direct measurement and by subse- 
quent analysis of some flows. 

An additional and very important influence on the flow in tidal basins is 
wind. Most catastrophic flood conditions in tidal basins that have resulted in 
levee breaks and other flood damage, were caused at least in part by wind. If 
a strong permanent wind flows over a water surface it creates a shear force 
which drives the water in the direction of the wind and makes it assume a 
steady slope towards the wind to counteract the shear force. Flow velocities 
will result as a consequence of changes of this wind force. Both the strength 
and direction of the wind can not be predicted with sufficient accuracy to in- 
troduce it as a cause of tidal calculations. Usually it is introduced only as a 
possible additional effect and taken care of by an additional stage and slope of 
the water surface, 

Other effects, such as temperature differences, which have an effect on the 
density of the water, may be assumed to be unimportant for the same reason 
as salinity differences as long as turbulent friction prevails. Animal and 
vegetable growth, on the other hand, may temporarily or permanently influ- 
ence both the friction factors and the location of sediment deposits and should 
be studied in important problems. 


The Differential Equations of Tidal Flow 


After defining a tidal flow in the preceding paragraphs and after discussing 
there the various effects and causes which must go into a mathematical de- 
scription of such flows, the actual equations can be established which describe 
the flow. Let us first see if there are any other problems which are solved 
and may lead the way for a general solution of the tidal flow problem. 

The “tidal flow” was defined as the flow in a basin, canal or canal system, 
which is at least at one point open to an ocean with known tidal action, which 
may or may not have other openings with steady or fluctuating discharges and 
which is often sufficiently large for its depth to develop significant friction. 
At the opening to the ocean the water surface elevation is varying in a peri- 
odic fashion due to tidal action. This will cause a wave action in the tidal 
basin originating from that point suggesting that the problem be treated as a 
wave problem. At the same time flows will occur between the ocean and the 
basin which will constantly change magnitude and direction and these flows 
will propagate into the basin—suggesting that the problem may be interpreted 
as one of flood routing, especially since the period of the tide is very long 
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compared with most other periodic waves. Fortunately the differential equa- 
tions describing these two problems are identical such that we are free to 
interpret the flow as a wave problem or as an unsteady flow problem. 

Before the differential equations of this unsteady or wave motion can be 
established it is necessary to classify the basin under consideration in a gen- 
eral manner. This classification is easiest understood if the problem is in- 
terpreted as a wave problem. Everybody has seen the wave pattern which is 
generated by a stone thrown into a still pond. Circular wave fronts will 
emerge from the point of origin and spread evenly over the entire water sur- 
face in all horizontal directions. Let us call such a pattern three-dimension- 
al since the wave, which is a displacement of the water surface in the vertical 
direction, spreads in both horizontal directions no matter how the horizontal 
axes are chosen. 

The two-dimensional wave motion on the other hand is one in which most 
of the wave propagation occurs in one preferred direction with all wave fronts 
normal to this direction. Such conditions are found predominantly in canals 
which are long compared to their width and where the source of the wave mo- 
tion is located at the end of the canal. Due to reflection at the side walls such 
waves will have the tendency to orient themselves such that the crests be- 
come straight and normal to the channel axis. This statement is only partial- 
ly correct for short waves, i.e. for wave lengths short with respect to the 
channel widths, but becomes more and more generally correct as the wave 
length becomes large with respect to the channel width. This statement calls 
for an estimate of the tidal wave lengths, for the purpose of comparison with 
the widths. The wave length may be derived as the product of wave period 
(something over 12 hours) and the wave velocity which for this purpose may 
be estimated by the shallow water wave propagation formula for frictionless 
flow, c =V@Y, where y is the still water depth and g the acceleration due to 
gravity. The wave lengths resulting for depths of 3, 10 and 50 feet are 85, 
150 and 340 miles, respectively. It will be agreed that these wave lengths are 
long with respect to the widths of most tidal basins; it will become rather in- 
teresting in the application of any resulting formulas that in most practical 
cases the wave length is of equal order of magnitude or even longer than the 
length of the tidal basin in which it propagates. The order of magnitude of 
wave period and wave length are the same for flash floods in medium sized 
streams and the resulting flows must be expected to cover also similar 
ranges of conditions. As all such flood flows are calculated by means of the 
hydraulic approach, it can be expected that the same approach, which essen- 
tially interprets the flow as two-dimensional, is also sufficient for many tidal 
problems. Up to this date only two-dimensional problems or those border- 
line cases which are three-dimensional but can be approximated with suffi- 
cient accuracy by a two-dimensional description have been studied analytical- 
ly in the U.S. 

In this section we thus derive equations for the velocity and surface eleva- 
tion of water in a channel as functions of position and time. These equations 
are not specifically concerned with tidal flow. The influence of the tide en- 
ters the solution of these equations through the boundary conditions imposed 
at the ocean entrances. It is clear that by introducing other boundary condi- 
tions one could in a similar manner solve problems such as flood routing. In 
fact some of the techniques of solution are applicable to quite general bound- 
ary conditions. It is only due to the approximate periodicity and the relative- 
ly small flow velocities of many tides that appreciable simplifications are 
possible, 
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In order to study the details of varied flow such as velocity distributions, 
one generally uses a hydrodynamic approach in which local values are solved 
for. In many cases a different approach is used either because specific de- 
tails are not important or because hydrodynamic solutions are not feasible. 
This latter approach we call hydraulic. The entire cross section is consid- 
ered as a stream tube of variable dimensions and the flow in it is described 
by mean values over a cross section only. One must be careful in carrying 
out the averaging, however. The best way of doing this is to proceed directly 
from the hydrodynamic equations and to integrate over a cross section. For 
such a derivation we refer the reader to a paper by Keulegan and Patterson. (2) 

We adopt a coordinate system in which the x-axis is horizontal, generally 
following the center line of the channel bed. The y and z axes are perpendic- 
ular to x and lie in the plane of the channel cross section, y being horizontal 
and z vertical. 

We shall assume, as in(2) that: 


1) The fluid is homogeneous. We specifically neglect thereby density cur- 
rents, sediment movement, etc. ‘ 

2) The fluid is incompressible. 

3) The channel slope i = tan @ is small so that tan@~ @= sin 9. 

4) The flow is gradually varied. This means that the second derivatives of 
the cross-section A and of the mean velocity U and the squares and 
products of their first derivatives may be neglected in comparison with 
the first derivatives. In general this requires that the banks of the 
channel be not much flatter than a 1:1 slope, that the slope of the wave 
profile << 1 and that the curvature of the wave profile < 2 Qu, a 8 
where u is the mean velocity over a cross section. In such cases we 
can neglect the components of acceleration parallel to a cross section. 


Let U (x,y,z,t) be the x-component of the local flow velocity and A (x,t) be 
the cross sectional area. We denote by 


Q= UdA = discharge 
A 
u =Q/A_ = mean velocity. 


Consider now the conservation of fluid for a small volume AV = A(x,t)Ax of 
cross sectional area A(x,t) and thickness Ax. (Fig. 1) 


z 


j-—— Ax 


Fig. 1: Channel Section. 


The net outflow of water in A V per unit time 
= Q(x+Ax)- Q(x) + ny 
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Since the fluid is incompressible, net outflow plus increase in volume equals 
zero. Thus 

a(uA) re) 

et (1) 


This equation applies to canals or estuaries alike if lateral flows are neg- 
lected. If, however, the local velocity changes very much over a cross sec- 
tion such as in an estuary, the quantity u loses its power as a parameter to 
describe the flow sufficiently. Of particular interest in such a connection is 
an estuary containing a narrow, deep channel and adjacent extensive shoals 
which are flooded when the tide rises. Another such case of interest is a 
river containing groins, the water circulating in eddies in the region between 
the groins in such a way as to guide the flow in the stream bed. In such cases 
it is convenient to divide the total cross section into a stream bed in which 
the flow velocity is relatively high and a lateral basin which contributes only 
to storage and not to discharge. In general the water velocities and water 
levels are quite different in the stream bed and the lateral basins. If we de- 
note by A the area of the stream bed and Ap the area of the lateral basin, the 
general continuity equation is (see Fig. 2). 


Assume that we are dealing with a situation such as that pictured below: 


Fig. 2: Cross Section 
If we assume that the water level has time to adjust laterally, then 


B _ dh 


A 
where the average width is bp = = Assume that b) is a constant so that the 
lateral basin has a rectangular cross section. Then 


_ ah 
at 


Thus the continuity equation becomes 
ax t dt ot 


where hy = b + bp. If moreover the stream bed has a rectangular cross 
section 115-7 


= 
The increase in volume per unit time 
aA 
Ax. 
bt 
b 
As' 
dhe ab 
at 


This is essentially the same as (1) except for the lateral discharge term and 
where instead of the width of the stream bed we use the total width by. 

Consider now the following contributions to the x-component of momentum 
per unit time for the volume AV: 


1) storage of momentum = Axp—*— f uda= 


2) net outflow of momentum = P GA = pax2-(aAu") 


u® 
uta 


3) force due to pressure gradient = —A xf dA=>A 


4) resistance force due to shear on wetted perimeter = -ax f tdx 
x 
where T = shearing force on wetted perimeter x. 


This can be written in the form Ax _ 
where C is defined by ro = -f — = 


We assume that the shear on the free surface is zero. 
We neglect any additional terms arising for estuaries with storage basins. 


For small velocities such terms are generally small enough to be neglected. 
Collecting results we have 


where R = A/x hydraulic radius. 
Neglecting the components of acceleration parallel to a cross section we 
find that the pressure p is given by the hydrostatic formula 


P = pg (H-z) (4) 
where H = surface elevation. This formula is valid for waves whose wave 
length is large compared to the depth of water (3). Substituting (4) into (3) we 
get 


et (uA) +3 (aAu )= Ag ax R (5) 


Multiplying the continuity equation (1) by ou and subtracting from the above 
we find 
ou 
+(I- - 6 


For a = 1 this reduces to 
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Sua), by 2h_-=0 (2) 

| 


du du _ gu lul 


The coefficient, a, appears to depend to a large extent on the width depth 
ratio. The numerical values are only slightly larger than one in most cases 
(see reference 4, p. 272). We have written u|u| in place of u* since for re- 
versing flows the frictional resistance must be proportional to u? and have the 
same sign as u. 

The hydraulic radius, R, is a function of x andt. To illustrate the varia- 
tion of R consider a canal of trapezoidal cross section as shown in Fig. 3. 


Fig. 3: Trapezoidal Section. 


COTa 
x 1+(2h/b)CSC a 


For an example we find for the Panama sea-level canal) that for tana = 1, 
b = 600 feet, h = 60 feet, the hydraulic radius is R = 0.86 h. Moreover, for 

a tidal range of 20 feet and water depths of 40, 50 and 60 feet, R is 36, 44 and 
52 feet respectively. 


Solution of the Differential Equations 


The differential equations (1) or (2) and (6) or (7) include most tidal and 
flood flows which can be described as hydraulic problems. There is no limi- 
tation on the shape of the channel other than its adaptability to hydraulic flow 
description and constant width over the range of tidal stages. There is no 
restriction on the shape of the tide or on the hydrographs which may create 
the flows and stages under consideration. Equations (1) or (2) and (6) or (7) 
are thus so general and all-inclusive that their knowledge is in itself not a 
great help. Only a solution of these equations which satisfies the boundary 
conditions of the particular channel and of the given tide cycle and hydro- 
graphs, gives actually the water surface lines and flow velocities which the 
engineer needs for the design of the canal and its works. To find the explicit 
solution of the equations for a given set of boundary conditions is in the gen- 
eral case impossible in closed form. Most generally applicable, but also 
most time consuming, are various numerical step-methods. To the writers’ 
knowledge, this approach has not been used in the U.S. to any significant 
amount in tidal flow calculations. 

The other method of approach is the introduction of simplifying assump- 
tions; the most common are: 


1) Average elevation of water surface constant along horizontal canal. 
2) Simple harmonic tide. 

3) Linear friction law. 

4) Constant rectangular channel cross-section. 
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5) Nonlinear terms are neglected in the differential equations. 


Various solutions are possible on this basis. Chronologically, the first 
was made by W.B. Parsons on te occasion of the solution of tidal flow prob- 
lems for the Cape Cod Canal. 6 


Parsons’ Harmonic Theory 


This is based on the introduction of all five approximations enumerated 
above. In this case an exact solution of the equations (1) and (7) is possible 
and was given for the case of two different tidal cycles applied to the two ends 
of a uniform horizontal canal. While approximations 1, 2, and 4 are easily 
understood, it may be necessary to explain the meaning of approximations 3 
and 5. It is clear that the turbulent friction in a canal will always be de- 
scribed by an energy slope which is about proportional to the square of the 
average velocity u. For a particular flow one may, however, approximate the 
term 


g lulu 
C*R 
where C' is a constant equal to a representative value of gluli/ C*R , 


such that the work done by the linearized friction equals that of the quadratic 
friction of equation (7) over a full cycle. This changes (7) into 
du du 


‘ 8 
at ay +9 (8) 


by C'u 


while equation (1) may be written as 
(9) 


since the width b was assumed to be constant over the range of tidal stages. 
For sinusoidal tides the principle of equivalent work gives 


3” 


The factor C' is in general a function of x, but is often averaged over a given 
reach for convenience of calculation. 

Both equations are non-linear in u, H and their derivatives. It can be 
shown, however, that these equations can be linearized for tidal flow problems 
if the small terms are eliminated. Such an appraisal can be based on the 
fact that all the waves under consideration are smooth and very long for their 
height. Let us introduce an auxiliary small constant parameter € which has 
the order of magnitude of the ratio of the wave height divided by the wave 
length of the tidal wave. 

Now let us develop the variables as follows 


U= + (10) 
H = H, + «eH, + «e*H, + 


Such Up and Ho, respectively, are the values of u and H for zero wave steep- 
ness. Introducing these expressions into (9), we find for equation (9) 
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Quo 4 
Uo Ste +H + 


This equation (11) must be fulfilled also for the case of zero wave, which is 
represented wr € = 0, and ren thus be divided into two independent equations 


(12) 
dH OHi 
Uo ° = 


In order to obtain we assume the average velocity 
in these linear equations between uj, Hj and their derivatives vanishes. 
With uo = 0, the first equation (12) gives dH,/dt = 0, which is correct any- 


way, since Ho is an average value, 
The first and the fourth term of the second equation disappear too, leaving 


of equation (12) 


aH ite aH 


Equation (8) may also be linearized by the use of (10) and gives 


(14) 
Ou, du Uo dH 


In the first equation (14) all terms except the third vanish because of 
U, = 0 - from which 3H,/dx = 0. In the second equation the terms with u,, 


disappear, leaving altogether from (13) and (14) 


if the indices | are dropped as insignificant. 
Ho is herein the constant undisturbed water depth. 

H_ is the deviation of the actual water surface from the zero level. 
u__ is the instantaneous water velocity. 


Equations (15) can be solved in closed form as Parsons showed, The solu- 
tion is rather easily obtained because of the linear character of equations (15) 
which permits the construction of the general solution by super-position of 
particular solutions. Parsons shows that the general solution for tidal flows 
in a horizontal canal with sinusoidal tides of equal frequency, but different 
phase and amplitude acting at the two ends may be given as the sum of four 
progressive damped waves. For the mathematical expressions, the reader 
may be referred to Parsons’ paper, since a duplication of the long expres- 
sions would serve no useful purpose at this place. 

Parsons applied his solution to the Cape Cod Canal and, after judicious 
choice of an appropriate friction factor C', was able to check rather closely 
the measured water profiles along the canal for various phases of the tidal 
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Heo + =O 
‘ (15) 
+9-$8 + cuz o 


motion. The proper choice of the friction factor is most important if a reli- 
able result is desired. A few remarks are thus in place about the choice of C'. 
The “linearized” friction factor C' is actually not a constant for any canal, 

but is changing in reality with the instantaneous flow velocity. If equation (15) 
applies an average of the instantaneous C' values, some of the details of the 
flow pattern may be lost, but its general shape is not much affected by this 
approach. 

The major advantage of equations (15) is their linear character which per- 
mits the superposition of individual solutions. Here arises the question of 
what velocity must be used to determine C' for the individual solutions. There 
is the tendency of applying for each individual partial wave its own average 
velocity for the determination of its C' value (see Brown’s reflected wave the- 
ory, which is discussed next). But that is wrong. A reasonably correct re- 
sult can be obtained only if the resultant velocity of all component waves is 
used to derive C'. This can most easily be seen if a wave is thought to be 
composed of two identical waves with equal phase and period as the original 
wave, but of half the amplitude. All other terms being truly linear, add from 
the two partial waves to exactly the same amount as determined for the orig- 
inal wave, all except the friction term in which both the velocity and C', which 
itself contains u, are reduced to half the original value. The friction term re- 
duces thus to 1/4 of the original value, and the sum of the two waves still has 
a friction term which is only half of the original value—proving this method 
to be wrong. It is easily seen, on the other hand, that the full friction term is 
obtained if the resultant velocity is used in the determination of C'. This way 
the inherently non-linear character of the friction term is preserved. This 
way there exists a mutual interference between the individual partial solu- 
tions. The proper C’ value can in most cases be obtained only by trial and 
error. 

Fig. 4 gives for example the solution of equations (15) to the model of a 
sea-level Panama Canal as described by Meyers and Schultz. It may be seen 
there that for instance the flow velocity at one end may be calculated with 
good accuracy from the vertical tide patterns at the two ends and a properly 
chosen C -value. 

Several such calculations were made in the past and instantaneous water 
surfaces were determined. For all canals which were short compared to the 
wave length as estimated previously, the profiles with high surface slope 
turned out to be practically rectilinear. This fact gives Pillsbury a valuable 
starting point in his calculation method. 

In conclusion, it may be said that Parsons’ method is correct and reliable 
where it applies, i.e. for sinusoidal tide applied to uniform horizontal canals 
without average flow in one direction. Whenever any of these conditions are 
significantly violated, one must change the form of the differential equations— 
and naturally of the solutions. 


Brown’s Reflected Wave Theory 


Tidal flow problems which are described by linear differential equations 
can be solved by superimposing the motion of simple progressive waves, and 
introducing the effect of abrupt changes in cross section. This is in essence 
the method proposed by Brown.(7) However, the actual solutions proposed by 
Brown appear unnecessarily complex and at the same time theoretically over- 
simplified. 


ae 
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u (ft/sec) 


TO PACIFIC — 


+ TO ATLANTIC 


@ MEASURED MODEL VELOCITIES 
(MEYERS AND SCHULTZ) 


CALCULATED VELOCTIES ACCORDING TO: 
—— PARSONS- FRICTION FACTOR FOR u = 3 ft/sec. 
—-— PARSONS-FRICTION FACTOR FOR u= 4ft/sec. 
—-— PILLSBURY 
—--—BROWN (ROSENZWEIG) 


Fig. 4 


Brown assumes that the effect of friction is small so that it can be consid- 
ered as a small change referred to a sinusoidal frictionless wave. The fric- 
tional resistance is taken as proportional to the product of particle velocity 
and wetted perimeter. Equating the work done by the frictional resistance to 
the loss of energy by a sinusoidal wave one finds, on integrating the resulting 
equation, a series approximation for the change of wave height with travel 
distance. 
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This derivation appears unnecessary, however, for assur_ing frictional 
resistance proportional to velocity the motion of a simple progressive wave 
can easily be determined exactly as demonstrated for instance by Parsons. 
The phase of the tide is there shown to be affected by the friction as well as 
the amplitude, thus contradicting Brown’s assumptions. Moreover, the flow 
velocity is not simply given in terms of the vertical tide by the formula for 
frictionless flow, as Brown supposes, but is given by more complicated 
expressions. 

Consider now the partial reflection of an incident progressive wave at a 
place of abrupt change in section of a canal from breadth and depth bj, Hj; to 
b2, Hg. The following determination of reflection at these changes of cross 
section is reasonable at distances from the change which are moderate mul- 
tiples of the transverse dimensions of the canal. Within this region we as- 
sume that the net effect of friction is negligible. Under these conditions the 
solution for the amplitudes of the transmitted and reflected waves is well 
known, (3, Sect. 176) The analysis given in (3) is essentially as follows. As- 
sume that at the change in cross-section mass is conserved and the wave sur- 
face is continuous. Denote by aj, ar, at the amplitudes of the incident, reflect- 
ed and transmitted waves respectively. For an incident wave moving from 
depth H; to H2 we have the well-known equations for a frictionless progres- 
sive tidal wave in a uniform horizontal channel with H the instantaneous local 
depth and Ho the average depth. 


H-H,= @ Sin(kx-wt) (16) 
and the corresponding velocity of the water particles 


u= 
On/ sin(kx—wt) = (H — He) 
and at x = 0 depth change actually occurs 


b,h, 
Ci 


where the transmission and reflection coefficients are respectively 


and 


These equations have the solution 

(19) 
One easily verifies that the sum of energies in the reflected and transmitted 
waves is equal to the energy in the incident wave, i.e. byjC1R? + bgC2T2 = 1. 


Although we have used periodic waves to compute R and T this is not essen- 
tial in the derivation (see ref. 3). 
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C, =/gh, » C, = ./gh, 
| 


As particular cases we see that when b2~0, R-+1 and when -1. 
Thus at a closed end a progressive wave reflects with the same sign and equal 
amplitude. At an open end it reflects with opposite sign and equal amplitude. 

Thus for a canal which consists of sections of constant cross section we 
have a complete scheme for dealing with linear wave motions. Within each 
individual section the roughness may vary. Starting with conditions at one end 
we proceed along the canal in small increments, advancing the solution step 
by step by Parsons’ equations for progressive waves. At a change in cross 
section, including an open or closed end we use formulae (19) to determine 
the new amplitudes and assume periods and phases are unaffected by the tran- 
sition, The resultant tides are simply formed by algebraic addition of suc- 
cessively reflected and transmitted waves due to tides acting at the two ends 
of the canal. If the canal has essentially constant characteristics explicit 
solutions are available as in Parsons. One should note that the friction coef- 
ficient must be estimated from the total velocity and not for each individual 
wave, as Brown proposes. 

Brown’s method is thus based on the very promising idea of treating the 
tidal flows as a system of progressive waves. The method has been applied 
to various existing canals and gave good results for the shape of the water 
surface profiles. It must be repeated that in most channels which are all 
short compared with the wave length, the profiles are usually straight lines. 
With the water surface elevations given at both ends, the profile is bound to be 
correct. Rather large deviations from measurement result in the velocities 
which depend heavily on the proper description of the frictional effects. But 
there is no reason why these shortcomings could not easily be eliminated. 


Pillsbury’s Theory 


Siase previously described methods appeared to be unsatisfactory, Pills- 
bury 8 proposed a method of successive approximations for dealing with tidal 
flows in canals and estuaries. The corrections are based on the so-called 
“primary currents.” 

The primary current is derived by first assuming that the vertical tide, H . 
has a known amplitude and phase while the flow velocity, u, has the same 
period as H‘) Thus assume that 


H'= H—Ho= h cos(ct+y) (20) 
u=B sin(ct + B) (21) 
where h,”, ¢ are known. Then by (20) the slope is of the form 


dH'_ 22 
ax = A cos(ct +a) (22) 


where A, ao can be computed from H, ¥. 
Substituting (21), (22) into (20) 


cos(ot +B) + gA cos(ct+a)+kot Bsin(ct+8)= O 
Equating coefficients of sin(gt) and cos(ot) separately to zero 


oBcos8+ gAcosa+k.$BsinB = O 
-7BsinB- gAsina+k, CB cosB= O 
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Solving these equations for B and 8 we find 


l-o cota 


(25) 


The variables A, @ are computed at the middle of the canal from the en- 
trance tides by assuming that the instantaneous profiles are straight lines. 
Thus denoting by Hj and Hj the entrance tides 


A cos(ot+a) (26) 


This assumption is nearly correct for short canals such as the Panama Canal 
(L~ 42 miles) as we noted previously. The amplitude and phase of the flow 
velocity at the middle of the canal are then computed from (24) (25). 

Flow velocities at other locations in the canal are obtained by summing 
known velocity increments. Thus writing the continuity equation in finite dif- 
ference form 

= Bk hg sin(ot+y) (27) 
We add these increments of velocity to the value calculated at the center of 
the canal in order to find the velocity at other locations. Thus this velocity 
will be of the form 


u=B sin(ot+f)= Bmsin(ot+Bm) +> Sine sin(ot+y)- (28) 
Hence 
B sinB = Bp sin he siny 
(29) 
Bcos 8 = Bm cosBm+5 ho cosy 


Further adjustments of the elevation differentials are made so that their sum 
is the difference of the entrance elevations. 

On the basis of the primary current computed in the manner described 
above, Pillsbury derives distortions due to the velocity head term in (7), the 
minor components of the friction and variations of the hydraulic radius, R, 
and Chezy coefficient K. Consider the correction i to the approximate veloc- 
ity B sin(@t +f) due to friction alone, i.e. let 


u*Bsin(ot+B)+i . (30) 
Neglecting the velocity head term in (7) and assuming a is given by (22) we 
find on substituting into (7) that 
(ot+8)+il= 
Now B, # are chosen so that the first equation of (29) is satisfied. Hence 
B + gAcos(ct+a) +k, B sin(ot+ 8) = (32) 
Thus by (31) and (32) 
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tan B = Lcota-o 
n(a—B) 


Expanding the last term in this equation, neglecting i* and also the effect of 
i on the sign of the velocity, one obtains a linear differential equation for i 
which is readily solved. Not realizing the possibility of an explicit solution, 
Pillsbury resorts to a laborious numerical solution by finite differences. In 
a similar way additional corrections can be obtained by setting 

u = B sin(ot + B) + i + j and expanding. 

This procedure of Pillsbury’s for computing corrections to a primary cur- 
rent independently for velocity, friction head, etc., does not appear correct. 
Moreover, and this seems to be the largest source of error, the elevation and 
slope can not be simply assumed to be of the form (20), (22) in computing 
these corrections. In fact, one of the principal effects of the nonlinear terms 
in the differential equations (2), (7) is to introduce additional harmonics in 
both H and u. 

Again, the same remarks can be made as for Brown’s method. The short- 
ness of most canals prevents any shortcomings of the method from influencing 
the water surface profiles much. Again, the velocities show the larger devia- 
tions, although an approximate choice of friction factors can go far in over- 
coming any basic difficulties of the method. 


Final Remarks 

It is remarkable how small the number of actually performed tidal flow 
calculations in the U.S. is, especially since work has been done on large num- 
bers of problems in this particular field. The writers have not been able to 
find any such calculations on estuaries, i.e. non-uniform channels. This may 
be explained, at least in part, by the large size of the projects involved: 
There existed a certain reluctance to base the large and important projects 
on approximate and partially oversimplified calculation methods as long as 
the solution by the “more trustworthy” model could be found. Th< writers do 
not want to close this paper without pointing to the main advantage of the cal- 
culation: even rather complicated calculations are many times less expen- 
sive than a model and can show many results at least as clear and distinct if 
executed by an experienced man or group. 
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